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I. INTRODUCTION 



Recently important progress has been made in the analytic treatment of self-gravitating spinning compact objects 
in general relativity using different methods. One of these methods is based on an effective field theory (EFT) point of 
view, which was applied, e.g., within the post-Newtonian (pN) approximation to non-spinning compact objects [1-5] 
^Jj- and to spinning objects [6-15]. An advantage of this approach is that some of the very sophisticated and systematic 
techniques for perturbative calculations used in high energy physics can be applied in a straightforward manner. The 
present paper addresses certain open questions regarding the handling of the spin supplementary condition (SSC) 
CO [ within this approach, though some aspects may be applicable to other approaches as well. In particular it is shown 
how the covariant SSC can be eliminated at the level of the potential (which is subtle in various respects) and how 
the dynamics can be cast into a fully reduced Hamiltonian form. For a review of spin in relativity and the problem 
^ : of the SSC see, e.g., [16-18]. 
^vq ■ In classical Newtonian mechanics the spin of an object is described by a 3-dimensional antisymmetric tensor or by 
— ^ , its dual vector. In some situations (e.g., when the dynamics depends on the spin but not on the absolute orientation 
I ■ of the objects) it is very convenient to associate a Poisson bracket representation of the so(3) Lie algebra (i.e., the 
] angular momentum algebra) with the spin variables and describe the spin dynamics in terms of a function generating 
,—1 . the time evolution via these Poisson brackets. This function may be a Hamiltonian or a Routhian [19, 20]. (The latter 
L| ' is a Hamiltonian for a part of the variables only and a Lagrangian for the remaining variables.) Such a formulation 
. ^ , of spin is also desirable in the relativistic case. But in this case the spin is given by a 4-dimensional antisymmetric 
■ tensor, so the best one can immediately achieve is to relate it to the so(l,3) Lie algebra of the S0(l,3) Lorentz 
, group. However, it is well-known that the fixation of a representative worldline or center of the spinning objects is 
- . ■ equivalent to a supplementary condition on the spin components. With this SSC the independent components of 
the 4-dimensional antisymmetric spin tensor are given by its 3-dimensional spatial part. But the reduction of the 
so(l,3) algebra for the spin components to a so(3) algebra is subtle and must be discussed within the framework 
of constrained Hamiltonian dynamics. Using the Dirac bracket approach this was performed for flat spacetimes [21] 
and for test-spinning objects in curved spacetime [22]. For self-gravitating spinning objects the reduction succeeded 
to linear order in spin with the help of an action principle [23, 24] and agrees with a construction via generators of 
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the global Poincare algebra valid to next-to-next-to- leading pN order [24-26]. But this derivation is focused on the 
canonical formalism of Arnowitt, Deser, and Misner (ADM) [27, 28] and is restricted to the Schwinger time gauge 
family of vierbein gauges [29] (which unfortunately does not include the gauge used in the EFT approach). 

A consistent way to deal with the (covariant) SSC within the EFT method is described in [8] (based on developments 
in [6], see also [13] and [7] for a sophisticated Kaluza-Klein-like reduction of the components of the metric tensor for 
further construction of the Feynman rules). There a Routhian generates the spin evolution via the so(l, 3) algebra and 
the center of mass motion is given by Euler-Lagrange equations (see also [30]). The covariant SSC is eliminated at the 
level of the equations of motion. In the present paper we show how the covariant SSC can consistently be eliminated 
already at the level of the potential using two different methods. One method is based on the well-known Dirac 
bracket and the other on an action principle. It is then straightforward to obtain a fully reduced Hamiltonian form 
of the dynamics. This is very convenient, e.g., for deriving the fully reduced equations of motion (where the length 
of the 3-dimensional spin vector is constant without further variable transformations) or for an implementation into 
the effective one body formalism, see [31-33] and references therein. We treat the pN next-to-leading order (NLO) 
spin-orbit, spin(l)-spin(2), and spin(l)-spin(l) level here. The relevant potentials were derived in [8-14] within the 
EFT formalism. Recently also the corresponding source multipole moments were calculated [34] (see also [35, 36]). 
The NLO spin-orbit and spin(l)-spin(2) dynamics was obtained earlier in [25, 37-40] and extended to arbitrary many 
objects in [41]. The NLO spin(l)-spin(l) Hamiltonian for binary black holes was given with correct center of mass 
motion in [42, 43]. In [44] the NLO spin(l)-spin(l) dynamics for general compact objects (including neutron stars) was 
reproduced and put into fully reduced Hamiltonian form. Higher orders in spin were treated in [43, 45]. For radiation- 
reaction effects on the motion of a binary due to spin see, e.g., [46] and references therein. Very recently even the 
next-to-next-to-leading order (NNLO) spin-orbit [47] and spin(l)-spin(2) level was tackled, the latter simultaneously 
by a potential within the EFT approach [15] and by a fully reduced Hamiltonian [48]. An extension of the results in 
the present article should be useful to relate these two results at NNLO spin(l)-spin(2). Notice that above results 
were obtained only very recently compared to the first treatments of self-gravitating spinning objects within the 
post-Minkowskian [49] and post-Newtonian [50, 51] approximations (see also [52-56]). 

It should be noted that the basic approach used in the present paper was already described in [44]. Therein also 
the result for the transformation to canonical variables was used in advance, but the presentation of the derivation 
was reserved for the present paper (some more details can also be found in [57]). In the meantime the basic approach 
was also applied in [9, 13], but still the transformation to canonical variables was not derived from general principles. 

The paper is organized as follows. In Sect. II an overview of the problem addressed in the present article is given. 
In Sect. HI the known potentials are transformed into non-reduced canonical form, i.e., still with a so(l,3) Poisson 
bracket algebra for the spins. In Sect. IV the Dirac bracket for the covariant SSC is calculated and transformed to 
standard canonical form, i.e., with a so(3) Poisson bracket algebra for the spins. An alternative elimination procedure 
via an action principle is performed in Sect. V. At the end of Sect. V it is discussed how this alternative approach 
can be used to improve the Feynman rules of the EFT formalism by formulating them in terms of reduced canonical 
spin variables. In Sect. VI the non-reduced Hamiltonians obtained earlier are transformed into fully reduced ones and 
compared with other results. Finally conclusions and outlook are given in Sect. VII. 

Our units are such that c — 1, where c is the velocity of light and also the implicit inverse pN expansion parameter 
with the formal counting rule 1/c" ~ fP-N order. Adapted to EFT convention we work in the spacetime signature 
—2 which is important to remember especially when working on the action level, see Sect. V. Three different frames 
are utilized in this article, denoted by different indices. Greek indices (a,/^, . . .) refer to the coordinate frame, lower 
case Latin indices from the beginning of the alphabet (a, b, . . .) belong to the local Lorentz frame, and upper case 
Latin indices from the beginning of the alphabet {A, B, . . .) denote the so called body-fixed Lorentz frame. Lower 
case Latin indices from the middle of the alphabet ■ ■ ■) are used for the spatial part of the mentioned frames and 
are running through (i — 1, 2, 3). In order to distinguish the three frames when splitting them into spatial and time 
part, we write a = (0), {i) for Lorentz indices (or a = (0), (1), (2), (3) in more detail), A = [0], [i] for the body-fixed 
frame, and /i = 0, i for the coordinate frame. Letters / and J are body labels, i.e. I,JG {1, 2} and z = (z*) denotes 
a point in the 3-dimensional Euclidean space M.^ endowed with a standard Euclidean metric (z* = Zi = S"^^ Zj) and 
a scalar product (denoted by a dot), so z/ e M.^ denotes the position of the Ith body. Indices appearing twice in 
a product are implicitly summed over its range, except for label indices of the objects. Round and square brackets 
are also used for index symmetrization and antisymmetrization, respectively, e.g., A^^'^^ = i (^'"^ + A'^'^). We also 
define r/ := z — z/, r/ := ]r/[, n/ := r//r/; and for I ^ J, r/j :— zj — zj, which is the distance vector of the two 
bodies. Likewise we define rjj :— [r/j], n/j :— r/j/r/j; ] • ] stands for the length of a vector. The linear momentum 
of the Ith body is denoted by p/ = (pn), and m/ denotes its mass parameter. An overdot, as in z, means the total 
time derivative. The spin vector of the Ith body is denoted by S/ and is always supposed to be referred to the local 
Lorentz frame in all potentials and Hamiltonians that are considered in the present article. The connection between 
the antisymmetric spin tensor S'(i)(j) and the spin vector is made by usage of the totally antisymmetric e-symbol 
^ijk = ^{i ^ " k){k — i) as Si^i)Q) = tijkS(k). Each body is within our approximation completely characterized by 
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its three parameters mass, momentum and spin. We associate to them the following relative pN order 

m.^o(l), p,^o(^l), Sj^O^^y (I.l) 

yielding the formal pN order for the Newtonian (N), leading-order (LO) and next-to- leading-order (NLO) Hamiltonians 
in question 

H^^O^^y Flo-o(^), ffNLO-o(^), (1.2) 

starting with the non-relativistic Newtonian Hamiltonian for two interacting bodies in canonical conjugate variables: 

p? pi Gmimz 
lm.1 2m2 ri2 

Notice that defines the zeroth pN order for Hamiltonians, so the pN orders of the LO and NLO Hamiltonians 
must be counted relative to iJjM (which is at 1 /c^ in the formal counting given above) . 



II. SETUP OF THE PROBLEM 



In [8] an EFT approach to include spin within the pN approximation is given, based on developments in [6] (see 
also [13]). This approach is able to deliver a pN approximate description of the dynamics of spinning compact objects 
via a Routhian, namely a function which is a Lagrangian for the objects 3-dimensional position z and a Hamiltonian 
for the 4-dimensional antisymmetric spin tensor Sab (in this section we may drop the object labels). Therefore this 
Routhian generates the spin evolution via an so(l, 3) Poisson bracket algebra and the center of mass motion is given by 
Euler-Lagrange equations (such an approach was already used in [30]). In order to allow a treatment of the covariant 
SSC via Dirac brackets one needs a canonical description of all variables, including z. Therefore we first perform a 
standard Legendre transformation in order to replace the velocity v = z by its generalized momentum p, which is 
of course the canonical conjugate of z. (Eventually accelerations or even higher time derivatives must be eliminated 
using the method in [58] first.) 

It is important that Sat is actually the generalized momentum of the 4-dimensional angular velocity tensor defined 

by 

^^°^ = AA°A'^^ (H.i) 

The latter is build from a Lorentz matrix A'^" G S0(1, 3), 

A^^A-^, = r^^b , (n.2) 

which relates the body-fixed and the local Lorentz frames; to be more precise A"*° is a representation of elements of 
S0(l,3). Here ijab = diag(l, — 1, — 1, — 1) is the Lorentz metric with signature —2. Notice that this property of the 
Lorentz matrix makes the angular velocity tensor antisymmetric. The complete Poisson brackets read, see, e.g., [21] 
or Sect. V below, 

{z\pj} ^ Sij , (II.3a) 
{A'^b,Scd}^VbcA^d-VbdA^c, (n.3b) 

{Sab, Scd} — ScaVbd — SdaVbc + SdbVac — ScbVad , (II. 3c) 

and all other zero. We will call these brackets canonical here, though they are not canonical in a strict sense (e.g., the 
spin components are not commuting). But it is possible to relate A"*^ and Scd to variables for which the phase space 
structure is manifest, see Sect. 3. A in [21]. For example, one can parametrize A"*" by independent angle-type variables. 
The time derivatives of these variables are angular velocities and their generalized momenta are canonically conjugate 
to the angle-type variables. (These angular velocities are also contained in 17°'', but with prefactors depending on the 
angle- type variables.) Notice that [8] uses a different sign convention for the spin part of the Poisson bracket shown 
above. Further the Hamiltonian does actually not depend on A^*^ in our case, so one can ignore A^° and its Poisson 
brackets for obtaining the equations of motion. For the present article, however, A"^° is of crucial importance. 

The Poisson brackets above are fully canonical, but the degrees of freedom are not fully reduced. An SSC must be 
imposed, which corresponds to a choice for the representative worldline of the compact object. A similar condition 
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must be given for A"^° [21]. According to [21] the covariant supplementary conditions read (consistent with the SSC 
used in [8]) 

5''V?=0, jyM^^A^f^, (II.4) 

mo 

where is Dixon's momentum of the compact object and rriu the dynamical mass. It holds m|j = g^'^p^p^ , where 
is the 4-diniensional metric. From (5.13) in [24] we get 

pf^ = + O {S^) , m = mD + OiS^), (II.5) 

where m is the constant mass parameter of the object. It should be noted that the mass-shell constraint is already 
implicitly eliminated within the approach in [8], as a gauge-fixing for the worldline parameter was performed. (Indeed, 
only a 3-dimensional canonical momentum pi is defined by the Legendre transformation mentioned above, but not its 
time component po-) The worldline parameter was chosen to be the coordinate time in [8], w° = 1. Notice that (II. 4) 
guarantees that in the rest frame the spin tensor contains the 3-dimensional spin S^-^^^^^ only (i.e., the mass dipole 
part S'(°)(*) vanishes) and that in the rest frame A'^^ describes a pure 3-dimensional rotation (no Lorentz boosts). 
This obviously reduces the degrees of freedom to the physically relevant ones, which are given by ^W^J) and AW(j). 

The most prominent way to handle the constraints (II. 4) on the phase space described by (11.3b, II. 3c) is provided 
by the Dirac bracket, denoted by {^-jo here. It is straightforward to calculate the Dirac bracket for the current 
situation in a pN approximate way. Notice that only three components of each condition in (II. 4) are independent, 
so one has six independent constraints for each particle. To the considered approximation the derivation will turn 
out to be very similar to [21], with the notable exception that in [21] the mass-shell constraint together with the 
gauge-fixing of the worldline parameter was also treated using the Dirac bracket. The Dirac bracket is essentially the 
Poisson brackets of the reduced phase space, but the variables used above are not canonical any more with respect to 
the Dirac bracket. Our next step is thus to transform z*, pi, S'(*)(-'), and AI*!'-'^ to new canonical variables denoted by 
a hat such that 

{z\p,}^^S), (II.6a) 
{AWO), ^ -5j,„aW(") -t- Jj^jAW^") , (II.6b) 

{^«(^), ^WW}d - (^.i^WC^) - + <5.fc5(^)(') - (5,,^(-'")W , (II.6c) 

and all other zero. The algebra for the spin was reduced from so(l, 3) to so(3) and the Lorentz matrix A"^" £ S0(1, 3) 
was transformed into AI*!'^''^ S S0(3). These variables are a suitable generalization of the Newton- Wigncr variables 
defined in flat spacetime [59], but they are not canonically equivalent to the Newton- Wigner variables introduced in 
[10] (for a discussion see [60]). Notice that the transformation to canonical variables is highly ambiguous, i.e., one 
may always perform a canonical transformation. It turns out that we are able to choose pi = pi here. Further we will 
actually not derive the transformation to A^^l^-'^ via Dirac brackets, as the Hamiltonian does not depend on it anyway. 
For the same reason the components of the spin tensor in the body-fixed frame S'jjj and thus the spin length s given 

by 2s2 = S'f.jfjjS'WIJl = S'(*)(J)5'(*)(J) are conserved. It should be noted that in [22] a canonical Newton- Wigner SSC 
for test-bodies in curved spacetime was handled by Dirac brackets directly (and consistently implemented into the 
action). However, from a general point of view it is very convenient to start with a covariant SSC, as this manifestly 
displays the covariance of the effective theory, in particular when higher dimensional operators are included in the 
worldline action. In the following the SSC is always assumed to be the covariant one, if not otherwise stated. 

Spin in relativity can also be treated by an action principle [61-63], see also [6, 21, 23, 24, 64] and appendix A of 
[65]. It is indeed possible to derive the transformation to reduced canonical variables using an action approach. The 
Poisson brackets (II. 3b, II. 3c) are essentially represented by a term in the action of the form piz'^ + ^Sab^^"" ■ After 
the supplementary conditions (II. 4) are inserted, one must find new variables such that this term takes on the form 

PiZ^ + \k\"''\ which precisely represents the reduced brackets (II. 6b, II. 6c). Here it is important that 

f^Wib) g S0(3) must be a 3-dimensional rotation matrix, A^'^l^*' A^'^l^-'^ = d^j. This approach is very similar to [23, 24]. 
Notice that one needs the transformation from A^*!'^-'^ to AI*!'^-') for the action approach, which is not necessary for the 
Dirac bracket approach. 
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III. LEGENDRE TRANSFORMATION 



The effective potential usually depends on velocities and positions when referred to a Lagrangian L defined as the 
difference between the non-relativistic Newtonian kinetic part Tn and the effective potential Vcs'- 

Leff = Tn - Feff = + - Kff . (III.l) 

The effective potential is the only part of the Lagrangian which is pN expanded and is further decomposed into 
different spin contributions. The first step to arrive at a reduced canonical Hamiltonian when starting with a non- 
reduced effective potential V^s = 14ff(x/, v/, 5'7(i)(j), 5'7(o)(i)) in pN approximation is to Legendre transform it only 
with respect to the velocities/momenta (the spin variables are formally kept unchanged by this procedure) to a non- 
reduced effective Hamiltonian i?eff = -ffeff(x/, p/, ^/(j)^), 5'/(o)(i)). In both expressions the SSC is not yet imposed, so 
that 5*7(0) (i) "will be treated as an independent variable. Furthermore the effective Hamiltonians are not supposed to 
contain any time derivatives of the variables except for the one of the position variable being defined as the velocity. 
If the case arose that a spin variable (including the constrained 5'7(o)(i) ones) would carry a time derivative one could 
either replace it with its lower order equations of motion or one could shift it onto positions and/or velocities in the 
same term by neglecting total time derivatives which serve as surface terms in the action. The last procedure ensures 
to leave us only with time derivatives of variables which are not further subject to a constraint, because the mass-shell 
constraint is already eliminated, when performing Legendre transformation. Those variables can therefore be treated 
differently aside from the Dirac bracket formalism but rather with a fully reduced Poisson bracket and a subleading 
potential/Hamiltonian, see the discussion at the end of Sect. V. But nevertheless one ends up with inserting lower 
order equations of motion for eliminating higher order time derivatives as outlined in [58, 66]. The Dirac bracket 
formalism will be carried out below for the SSCs in order to find a canonical set of variables after using the SSCs. 
The effective potential VcS for two interacting bodies is pN expanded up to NLO spin effects: 

Vpp is the point particle interaction potential. This again is decomposed into 

Vpp = Vn + VeiH/IpN + ^2pN , (III.3) 



starting with the Newtonian potential 



F. = -^™, (III.4) 



ri2 

and continuing with the Einstein- Infeld- Hoffmann (EIH) potential (Veih/IpN = — ^eih) [1, 67] 

1 ^ — T 4 GT]fl\lfin2 / 9 9\ / \ / \ / \"i 

^EIH = o / '"^a.-^a ^ 7, Vl + ^2) - 7 (vi • V2) - (vi • 1112) (V2 • 1112)] 

O — ^''12 

a 

G^mim2{mi + 1712) 
^'12 



(IH.5) 



The 2pN point particle potential is in fact not needed for the Legendre transformation, the interested reader can 
find it in [2]. The leading order (LO) spin contributions are decomposed into the LO spin-orbit (SO) contribution in 
non-reduced form, see [12]: 



LO _ Gm2 , 

^SO - 9 "12 
'12 



^W(o) ^ ^p)W _ 2«2^) +(lo2), (III.6) 



the LO spin(l)-spin(2) contribution, e.g. from 



= 3- [3(Si • ni2)(S2 • ni2) - (Si • S2)] (IIL7) 



' 12 

and the LO spin(I)-spin(I) (finite size) contributions from [10, 36] 
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with the spin quadrupole constant Cq^ which is chosen such that it is equal to one for black holes and correspondingly 
bigger for white dwarfs or neutron stars. The Legendre transformation with respect to the velocities/momenta is done 
by using the formula 



-ffcff = Vi • Pl + V2 • P2 - icff 

= vi • Pl + V2 • P2 - ^"iiVi - ^m2vl + VcS 



(IIL9) 



Pi 



1 



2mi 2m2 2mi 



(pi - mivi) 



(P2 - m-2V2)^ 



2m2 



with 



Pi 



dvj 



dV, 



This leaves us with the expression 





Pi 


' i 


2mi 


2m2 


2mi \ 




P2 ^ 


H{dv) 


2mi 


2m2 


1 




2 



off 



9v/ 



1 fdV,s 



2to2 V ^'^2 



Kff 



1 



2mi \ dvi 



dV,s 



2m2 \ dv' 



(III.IO) 

(III.ll) 
(III.12) 
(III.13) 



with V/ — V/ (p/,p,7)ipNj so it is sufficient to know the momentum up to IpN order, because at 2pN level correction 
terms to the Hamiltonian will be induced only by the square of the velocity derivative of the potential. The momentum 
is given by 



Pi = mivi 



d{V, 



EIH/lpN 



9v/ 



1 + -Vj ) TO/V/ 



Gmimj 
2rij 



[6v/ - 7vj - (n/j • vj)n7,j 



G 

+ —[mj{nij X Si) + 2m/(n/j x Sj)] + 0{c-^) . 
'''ij 



(III. 14a) 
(III.14b) 



The derivative of the potential is given by 



off 



9v/ 



nil 2 



Gmimj 
2rij 



[6v/ - 7v,7 - (n/j • v,7)n7j] 



G 



[mj{nij X Si) + 2TO/(n/j x Sj)] 



' ij 



(III.15) 



+ 0{c-') 



After evaluating its square and replacing the velocity by inverting Eq. (III. 14b) 



'Pi yPJ (n/j-pj) 
6 7 nij 



H 


4) 


P. 


Gmj 


mjj 


mi 


2rij 



mi 



mj 



mj 



G 



• IJ 



mi 



(n/j X Si) + 2(n/j x Sj) 



(III.16) 



+ Oic-'), 
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the Legendre transformation can be performed. The contribution Hf^gyy from (III. 13) then reads 

G ( 1 



^ ni2 • (pi X Si) r^ni2 • (P2 X Si) 



7TOini2 • (p2 X S2) + 6m2ni2 • (pi x S2) 

S^ + (Si • ni2)' ) + 2mi f - S^ + (S2 • ni2)^ 



(III.17) 



' 12 



2mi 



2m2 -(Si-S2) + (Si-ni2)(S2-ni2) 



+ (lo2)+0(c-i°). 



Certainly the replacement of the velocities through (III. 16) has also to be done in the LO spin-orbit and the EIH 
potential to arrive at the fully correct NLO Hamiltonians. The LO spin-orbit potential (III. 6) yields the contribution 
H^QQQ excluding the SSC, which remains still untouched: 



'loso 



G 




— 2 


Si 










''12 




- m2 





m2 
mi 



TO2 7 -I- 



_pL 

2mj 



(ni2 X pi) - 1 2 (ni2 X p2 



3m2 
mi 



Si • (ni2 X pi) + (6mi + 7m2) Si • (1112 x P2) 



m2 
mi 



(III.18) 



(Si -82) 



Likewise the EIH potential gives rise to the contribution H^'^- 



ttV^P 
^EIH 



G I 1 2// i~i \ \ ^2 2 

,p2((p2xSi).ni2)-^p2 
3m2 



' 12 
G2 



-7m2 



mi 

+ (lo2) + 0(c-i°). 



2ml 

2 ' ((pi X Si) • ni2 



pf((pi X Si) -1112) 
7m2\ 



6mi 



((P2 X Si) • ni2) 



(III.19) 



Now we are able to evaluate from the effective potentials for the NLO spin-orbit, spin(l)-spin(2) and spin(l)-spin(l) 
case their effective Hamiltonian counterparts Hcs — -ffcff(x/, p/, S'/(i)(j), iS'7(o)(i)), which result from Legendre trans- 
formation with respect to the velocities/momenta only. So these Hamiltonians still remain non-reduced in phase space 
as long as the SSC is not imposed. 



A. The non-reduced NLO spin-orbit Hamiltonians 



We consider two effective potentials from literature, one from Levi [13], and the other one from Porto [12]. Both will 
be subject to a Legendre transformation to arrive at -ffeff(x7, p/, 5'/(i)(j) , 5'/(o)(i)). These Hamiltonians still depend 
on the SSC, so they are far from being canonical and deserve only formally to be called Hamiltonians, in the sense 
that they are the result of a Legendre transformation of the potentials but only to a subset of variables, likewise the 
reduction of unphysical degrees of freedom is only achieved on a subspace of phase space. Both Levi and Porto include 
in their NLO potentials the SSC term arising from the LO SO potential, because the elimination of the SSC is itself 
subject to a pN expanded expression and will therefore lift LO expressions in the potentials to NLO ones and so one. 
For this reason it is important to keep track of the SSC in all terms where it is present. 



8 



1. The non-reduced NLO spin-orbit Hamiltonian of the potential of Levi 

We Start with the effective NLO SO potential of Levi as given in [13] Eq. (109) {V^q^) = --^so(l)): 

'1 1,3 



rNLO _ Gm2 ^ 
^SO(L) - — 2— »1 
'12 



Vl X ni2 ( |vi • V2 - ^vl - |(vi • ni2)(v2 • 1112' 



+ V2 X ni2 (vi • V2 - W2 + 3(vi • ni2)(v2 • ni2)) + Vi x V2 -Vl • ni2 + V2 • ni 



Si 



' 12 



Vl X ni2 ( 2toi - -m2 ) + V2 x ni2(2m2) 



(in.20a) 



Gm2 



Sf^^'^nl^ ( 1 - ^vi • V2 + "^vl - ^(vi • ni2)(v2 • 1112) 





^ Gto2 




(-^vi.ni2) 






ri2 


2 ^ ^ 



€"^7712 ^(0)(^) 



^I'^^^'^nl^ [mi + 2m2] 



' 12 



This potential owns a term with a time derivative of ^i*^^ According to our agreed rule we shift it onto positions 
and velocities in the same term yielding 



rNLO _ G7ri2 
^SO(L) - — 2~^1 
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Vl X ni2 ( -vi • V2 - -vl - -(vi • ni2)(v2 • 1112' 



V2 X ni2 (vi • V2 - W2 + 3(vi • ni2)(v2 • ni2)) + vi x V2 ^^vi • ni2 + V2 • ni2^ 



G^m2 



r' 



Si 



12 

(0)(i), 



vi X ni2 ( 2?7ii - -m2 ) + V2 x 1112(27712) 



' 12 



"^1 ^2 1 ^2 ' 



(in.20b) 



+ ^r^^'^77l2 ( 1 - ^Vi . V2 + ^7^2' - ^(Vl • ni2)(v2 • ni2) 



^G7Ti2S'i' 



(o)W 



a|_ _ (vi ■ ni2)7;2 (v2 ■ ni2)7;^ 
ri2 



12 



' 12 



G27772 



-^^f ^<*^77l2 [mi + 27^2] 
''12 



Next we insert lower order equations of motion to eliminate the accelaration term. By counting the pN order of this 
term in the potential it is clear that only the Newtonian equation of motion is needed as replacement which reads 



a2 



G777ini2 



(in.2i) 



'12 



The resulting potential will then be subject to a Legendre transformation meaning we replace velocities in (IIL20b) 
by momenta indicated by {v — ?> p) in ^g'o^(L)''^^'' and add the specific NLO spin-orbit contributions from (IIL17), 
(in. 18) and (III. 19), indicated by ~, to arrive at the effective Hamiltonian ff^Qj^^^*''^'' (x/, p/, 5/afc, 5'/(o)(i)): 

(III.22) 



„NLO(eff) „ I fT-i'^P I fT-f^P I -J .NLO(i;^p) 

■"SO(L) - tl(dVY + ^LOSO + ^EIH + ^SO{l.) 



leading to the result 



NLO(cff) _ 
SO(L) ~ 2 



^^P?((Pi X Si) • ni2) - ^^(pi-ni2)(p2-ni2)((pi x Si) • riis) 



+ ^(Pi-P2)((pi X Si) • ni2) 



1 



2toiTO2 



P2((Pi X Si) • ni2) 



3 1 

{Pi-ni2)(p2-ni2)((p2 X Si) • ni2) H (Pi-P2)((P2 x Si) • 1112) 

mim2 miTO2 



1 



2ml 



(Pi-ni2)((pi X Si) • P2) 



1 



mim2 



-(p2-ni2)((pi X Si) • P2) 



^( -5m2((pixSi).„i2) 



^ ((pi X Si) • ni2) + 6mi((p2 x Si) • 1112) 



llm2 \ G 

^((P2XS1).„12)1-^ 



-^1 "12 I - ' 2 2 



^(pi . ni2)(p2 • ni2)] + 5i°)(S2 (-^Pi • ni2 



+ ^^S[°^^'^n\^ [mi + 2m2] 



' 12 



Gmim2n\2 (pi • ni2)p| (p2 • ni2)p2 



' 12 



miri2 



m2r^2 



(III.23) 



^. r/ie non-reduced NLO sptn-orbtt Hamtltonian of the potential of Porto 
The second alternative potential we find in Porto [12] Eq. (53) 



Gm2 



NLO 

SO(P) - 2 
'12 



:W(o) /■ 1 , o,,2 o,, „ 3^^^ ^2 G , 



S\'""' \ l + 2v^ - 2vi • V2 - -(V2 • ni2)" - -(3mi + 2m2) 

1 - ^(V2 • ni2)2 + g (4mi - m2)) S^^^Vj 

2 - 2vi • V2 - 3(v2 • ni2)2 + ^v^ - ^ (2mi + 5m2)') ^^■' VAn 



+ Di [Wi — V2j V2 • ni2 + Di Vj^V2V2 • ni2 



+ (1^2) 



(III.24) 
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Notice that accelerations were already eliminated by inserting equations of motion. Legendre transformation yields 
like in the case of Levi's potential: 



TNLO(off) 
^SO(P) 



iier^.D^ ^ + -HlOSO + ^EIH + ^SO{F) 



NLO(eff) 
SO(P) 

3 



G 



' 12 



v—^p j_ -(^NLO(d— !-p) 

3 



Pi((Pi X Si) • ni2) 



2mf 



2toiTO2 



(p2-ni2)^((pi X Si) • ni2) 



rrin 



(p2-ni2) ((P2 X Si) • ni2) 



mim2 



(Pl-P2)((P2 X Si) • ni2) 



^P2((P2 X Si) • ni2) + — - — (p2-ni2)((pi X Si) • P2) 

1712 miTO2 

-m2 + ^(P2-ni2)^ + :^(Pi-P2) - ^pI] 



mi 



7712 



12 



— (,P2-ni2jOi Pi H iP2-ni2jO;^ P2 

7771 



G2 



57772((P1 X Si) • ni2) 



7777^ 

2-((pi X Si) • ni2) + 7777i((p2 X Si) • ni2) 



27771 



+ 67772((P2 X Si) • ni2) + 3777l77725r^^'^77l2 + 2777^ 5}"^ ^'^Tll 



(III.25) 



B. The non-reduced NLO spin(l)-spin(2) Hamiltonian 

We take the NLO spin(l)-spin(2) potential V^5°p) of Porto/Rothstein from [8] Eq. (56) (modulo the non-SSC- 
dependent LO spin-orbit terms while keeping the important SSC dependent LO spin-orbit term) 



NLO 



SlS2(P) 



G_ 

' 12 



{S'^ - 37712^2) ( S['^^°^ Si'^^°^ + ^Vi • W2Sr''''S., 



2 



-I- 3vi • ni2(w2r42 + ■f^2'^i2) + 3v2 ' ni2{v{n\2 + vlnl^) - v\vl - v\v{ 



(377'i2V2 • ni2 - 7;^)5(°)('=^5f + [Zn^^^i ■ ni2 - v[)S'^'"^^> S\ 

\^12 ^12 / ^12 



iAc(0)(/c) Mil) 



J (j)(0) 
„2 "'12"^2 
'12 



(in.26) 
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and apply to it the same Legendre transformation procedure as in the preceding section yielding 



^NLO(off) ^ jj 



'SlS2(P) 

rrNLO(cff) _ 

-"SlS2(P) ~ 



G f 15 

r ^ Tr(Pi'"i2)(P2-ni2)(Si-ni2)(S2-ni2) 

mim2ri2 \ ^ 



3 9 

2(Pi-P2)(Si-ni2)(S2-ni2) + -(p2-ni2)(Si-pi)(S2-ni2) 

3 3 3 

-(p2-ni2)(Si-ni2)(S2-pi) - -(pi-ni2)(Si-p2)(S2-ni2) + -(Si •p2)(S2 -pi) 

9 5 
-(pi-ni2)(Si-ni2)(S2-p2) - -(Si-pi)(S2-p2) 

3 3 \ 

-(pi-ni2)(p2-ni2)(Si-S2) + -(pi-p2)(Si-S2) I 



G 



' 12 



((Pi X S2) • ni2) ((p2 X S2) • ni2) S{ 



nil 



1712 



'12 



— {P2-ni2)eijkn\2S!f^ 3^°^^''^ —ei-jkPiiS2(,)S[°^^''^ + —eijkP2iS2U)S^ 



— ((P2 X Si) • ni2) - — ((pi X Si) • ni2) ) S^"""n\2 



(0)(fc) 



(0)(«)^-* Q(0)(j)^i 



OOi "■i2"-'2 

3 



.(o)W c.(o)W 



12 '^l '-'2 



— (,Pi-ni2jei77c«i2'Ji '32 H eiifePi»Ji '-'2 

mi -L ^ J ± ^ 



'^"^2 , <.0)(0) _ Gmi ■ (j)(o) 

2 12 1 li^'yiD'y 



'12 '12 

+ — 7(mi +m2)(Si-S2) - 13(mi + m2)(Si •ni2)(S2 -1112" 
''12 \ 

where ~ indicates here focusing only on NLO spin(l)-spin(2) terms with the SSC untouched. 



(III.27) 



C. The non-reduced NLO spin(l)-spin(l) Hamiltonian 

We adopt the NLO spin(l)-spin(l) potential of Porto and Rothstein that was calculated in [10] up to a missing 
contribution stemming from an acceleration term, which was corrected in [11]. The potential reads according to Eq. 
(49) of the arXiv version with the LO spin-orbit SSC term included: 



NLO 



CQi 



Gm2 



^?(P)"^^^2mir3 
- 25f ((vi X Si)'^ - 3(n • vi)(n x Si)'^) 



Gm2 



Sl [ 6(n • wif - yn • vin • va + y vi • V2 - - ^vj - 2ai • r 



(Si • nf 0(v^ + wl) - y vi • V2 - yn • vin • V2^ + 2vi • SiVi • Si 
3vi • S1V2 • Si — 6n • vin • SiVi • Si + 9n • van • SiVi • Si 



+ 3n • vin • S1V2 • Si 

G2m2 



Co, „ , i H 



TOl 



(S? - 3(Si . nf) 



(Si • n) V fa^i^i)) 5^°^('^+vixSi-a' 



Gm-, 



1(1) ^ ^2 
'12 



2„J q(j)(0) 



'n2'3i 



(III.28) 



12 



with the Si-dependent part of the acceleration of the local frame 



1(1) ~ T 



[-3v X Si + 6n(v X Si) • n + 3n • v(n x Si)] 



The acceleration term with ai appearing is eliminated by using the Newtonian EOMs for two bodies: 

Gm2ni2 
ai = 5 . 



(III. 29) 



(III.30) 



This potential will be Legendre transformed like the other potentials above resulting in the effective SSC-dependent 
Hamiltonian (~ indicates here sole focus on spin(l)-spin(l) terms with the SSC-dependent terms untouched) 



„NLO(off) rr , ttV^p tt^^P , taNLO(d-s-p) 

Sl(V) - ^(dV)^ + ^LOSO + ^EIH + ^ Sl(V) 



which results in 

fNLO(eff) _ G'^m2 



H 



S?(P) 



+ ( 2 + 3Gq, 



m2 
mi 



3Gr 



TO2 



o+6Cq, — (Si-n 



mi 



112; 



_G 

^12 



^((Pi X Si) . n,2)SrH. - -((P2 X Si) . ni2)5r«nl2 + f^l' 
mi mi 2mi V 



Cq, m2 _g(o)Mg(o)M 
2 mi 1 1 

(o)W^» , 3GQ,m2 /c,(o)(j)^ 



'12 



3(1 - GQjm2 



(pi-ni2) (p2-ni2) eyfeni2S'i(j)5i 



mi 



:(0)(fc) 



0)c.(0)(fc) 



„2 "'12'^1 
'12 



+ ( f-6 + ) P?(Si-ni2)' + (9- 3GqJ (pi •ni2)(Si •ni2)(Si -pi) 



(III.31) 



(III.32) 



+ (-3 + Gq J (Si -pi)^ + (-3 + 3Gq J (pi •ni2)^Sf + 3 



7Gr 



+ 



G 



'1' 12 



^Gqi (pi •ni2)(p2 •ni2)(Si •ni2)^ + f 6 + ^Gq, ) (pi •p2)(Si -1112) 



3- 2^Qij (P2-ni2)(Si-ni2)(Si-pi) - - (Pi •ni2)(Si •ni2)(Si •p2 

3-^Gq,) (Si.pi)(Si.p2)+ ("3- ^Gq,) (pi.ni2)(p2-ni2)S? 



Qi „2 



-3+^GQ,)(pi.p2)Sn+^^. 
4 ^ / / mim2rj2 



2 3 2 



p^ -(Si.ni2)^--St 



Notice that while the LO potential (III. 8) is purely spin quadrupole dependent (via the constant Cq^ ) the correspond- 
ing NLO potential (III. 28) is not, likewise in the case of the Hamiltonian (III. 32). 

These non-reduced Hamiltonians are now fit for further phase-space reduction procedures, either by Dirac brackets 
or by the action principle in Eq. (V.l). 



IV. REDUCTION VIA DIRAC BRACKETS 



Canonical formalisms in the presence of constraints were analyzed in a very general way by Dirac [68-71], for 
reviews see also [21, 72-74]. A very important tool developed in this area is nowadays called the Dirac bracket 
and is further explained in the following. Other important contributions to constrained Hamiltonian dynamics were 
made by Bergmann and his collaborators, e.g., the notion of primary constraints and the understanding of gauge 
transformations [75]. Early contributions were already made by Rosenfeld [76], e.g., the discovery of what is nowadays 
called the Dirac or total Hamiltonian. For a historical review see [77]. 



13 



A. Construction of the Dirac bracket 



The Hamiltonian formulation of a dynamical system needs the construction of a Poisson bracket type structure of 
the dynamics. In the case of a constraint dynamical systems this is not an easy task, but the needed formalism is 
available as developed by Bergmann and Dirac. In the following we shall present the formalism to the extent we will 
need it. At the beginning let us treat the following variables (^z},pjj, S"^, ^l^) as unconstrained. The time derivative 
of any function Q of these variables reads 

If one wants to restrict the time derivative to the independent degrees of freedom in all variables, one still has to 
reduce to the six degrees of freedom in the Lorentz matrices, which can be achieved by writing (IV. 1) in the following 
way 



Taking into account, see Sect. V, 



dp 



li 



P^' = ~~d^' ^Iab^'iSj^laVb]d-g^^~2f]^la^f^-^^ , (IV.3) 



^ab \.\Bb_xB nab _ r,dHcS 



A^"=A^fc^]"^ AAbA^'^S^, (IV.4) 

where, as seen from the previous equations, generally Hcs = Hcs{zi ,pj , S"^ , Af^^) will hold, we get: 

Q = {Q,H,s}. (IV.5) 

The Poisson bracket {■, ■} as defined here has its standard properties which comprise bilinearity, fulfillment of the 
Leibniz rule and of the Jacobi identity when performed with standard canonical variables having vanishing Poisson 
brackets among themselves. Hereof one can derive a chain rule analogon for the Poisson bracket when applied to a 
continuous function H depending on a canonical variable ^ which reads 

{Q,i/(e)} = {Q,a^. (IV.6) 

Applying this formula one can read off all the Poisson brackets, see Sect. V for a thorough derivation of them. The 
results are (V.10)-(V.12) and read 

{z],PJJ}^S^jSIJ, (IV.7a) 

{Slab, Sled} = SicaVbd — SldaVbc + SjdbVac ~ SjcbVad , (IV. 7b) 

{Slab, Af^} = TJbcAf^ - VacAft, , (IV.7c) 

all other zero. The physical evolution of the variables needs a SSC. We use the covariant one by Tulczyjew written in 
the local frame, Eq. (II. 4), as this SSC is the one implemented in the potentials shown in the last section (i.e., this SSC 
is conserved by the time evolution given by the potentials). Referring to (II. 4) we may replace Dixon's 4-momentum 
with the 4- velocity in the SSCs due to their equivalence to our approximation, cf. (II. 5), while avoiding unnecessary 
confusion when dealing with Dixon's 4-momentum whose 3-components pf ~ J^^u-^ essentially different from the 

canonical ones defined by Legendre transformation pf = ^ . This means that for the considered approximation the 
Tulczyjew SSC is equivalent to the Mathisson-Pirani one, SabU^ = 0. With the 4-velocity we get 



Sabu'' = ^ %)f,w'' = ^ 5(,)(o) +S'wo)-(^ =0, (IV.l 



, U( 



^ aW(°) +AW(^')-^ =0, (IV.8b) 

U(0) 
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Let us call the SSCs, including the A-relations, <I>^ = 0,{A = 1, 2, 12). Then the Poisson brackets 

{<S>a,<^>b}^Cab (IV.9) 

with the inverse matrix C^^ , i.e. C^^Cdb = 5^, are most important. Notice, the assumed non-degeneracy of Cab 
even for (f>^ — makes the SSCs to be second class constraints. The Dirac brackets are defined in the form 

{Qi, Q2}d = {Qi, O2} - {Qi, ^aIC-^'^I^s, Q2} . (IV.IO) 
The constrained evolution for the independent variables reads 

^/ = {^1; Hcs}t> , Pli = {PU, ^^cffjo , = {Sl(i)(j),Hcs}T> ■ (IV.ll) 

The Dirac bracket therefore satisfies all the laws known from the Poisson bracket, which turns it also into a Lie bracket 
and it leads to the correct equations of motion together with the Hamiltonian -ffeff- The Dirac bracket can thus be used 
as substitute for the Poisson bracket. However, whereas one may use the constraints only after all Poisson brackets 
were calculated, the second class constraints ^a — can be used before an application of the Dirac bracket without 
changing the result, e.g. one has {Q, <I>q}d — for all Q's and <i>Q thus preserving the constraints in time when using 
the Hamiltonian for Q. Notice if one restricts to use the Dirac bracket instead of the Poisson bracket, the second 
class constraints $a = can be used off-shell to solve for certain phase space variables and eliminate them from all 
quantities, thus reducing the actual degrees of freedom. The transition to new variables S'^''^"''') which fulfill 

{z\,p.j,}^^5]5u, (IV.12a) 

S^lfU)^Sbn)(n)^^ ^ -(5,>„AW(") + Sjnkf^"^^ , (IV.12b) 

{5f 5f = S,,sf^^'^ + - SaS[^^^'^ , (IV.12C) 

all other brackets being zero, results in a standard canonical representation. They will be called Newton- Wigner 
variables, because it can be shown that the ones proposed by those two in [59] represent the only possible standard 
canonical set of variables at least in Special Relativity. For an extension of definition of those variables to General 
Relativity see the comment [60]. The 'unhatted' variables corresponding to the covariant SSC (IV.8a)-(rV.8c) are from 
now on shortly dubbed 'covariant' variables. More details are given in the application later on. The careful reader 
should have noticed that the space and momentum variables {z}, pjj) were assumed non-constrained throughout. This 
was done for simplicity reasons and because the pure spinless case is well known. On the other side, the presented 
form is just the one delivered by researchers using the EFT approach, see later on. 

Now the supplementary conditions (IV. 8a)- (IV. 8c) need to be explicitly known in dependence of positions and 
canonical momenta in coordinate space when eliminating them in the Hamiltonians. Upon agreement the spin is 
always defined in the local Lorentz frame, so what we need are the 4-velocities in the coordinate frame mapped into 
the local frame, which is linked to the coordinate one by a vierbein transformation matrix. By decomposing the 
curved spacetime metric according to 

5/^1/ = Vf^iy + h^i, (IV. 13) 

into the Minkowski background field part r/^^ and a perturbation part /i^i^, whose indices are raised and lowered with 
ry^jy, we can derive for m° the expanded expression 

< = e''^(z/X = 77«^ (^77^, + ^h^, - ^h^pirK, + ..^u'i (IV.14) 

up to a certain pN order. (Notice that a vierbein is not uniquely fixed by the metric, but the one used here is entering 
the derivation of the Feynman rules.) Expressed in variables defined in an Euclidean fiat space coordinate frame the 
vierbein components are given up to all terms necessary for our declared approximation according to [60] by (see also 
[78]) 

6(0)0 = 1 + J/loo = 1 - G— - G— + 0{c~^) , (IV.15a) 
6(0), = 6(.)o = W = 2G«1^ - G5«(^)4 + (1^2) + 0{c-^) , (IV.15b) 
Hi)o = -'^y + \h^o = -'^y (1 + + g'^\ + 0{c-^) . (IV.lSc) 
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The velocities depending on the momenta follow from Legendre transformation using the pN expanded Lagrangian 
for two spinning compact body interaction leading to Eq. (III. 16). Notice the vierbein is chosen to be symmetric so 
that local Lorentz indices and coordinate indices get indistinguishable yielding the 4- velocities in Euclidean flat space 
coordinates and depending on canonical momenta which replace the coordinate velocities via Eq. (III. 16) indicated 
by:.. 



,(0) 



(zi)< = e(o)o(zi)ii; + e(o).(zi)< = 1 - + 0{c-^) 



mi 



(ziX 

1 - 



P? 



o 2 

2mj 



e(j)o(zi)wi - e(.,)j(zi)uj , 

2Gm2Pii 3Gp2i , 6*^2 (P2 -1112) 



Gm2nf2S{ 



{i)(k) 



TOiri2 



2ri2 

+ 0(c-5). 



2ri2 



(IV.16a) 
(IV. 16b) 

(IV.16c) 



The covariant SSC (IV. 8a) is then pN expanded to yield 



= 5i 



(o)(») 



S\ 



mi) 



' 12 



,(0) 



Pil 
mi 



ri2 



^ qjk ^12 
'12 



2ml 



G 



2ri2 



s: 



m2g(o){k)\ 
mi ^ 



0{c-^) 



3p2j - 

0(c-«). 



m2 , 

2 Pij + (ni2 • P2)ni2 

mi 



(IV.17a) 
(IV. 17b) 



(IV.lTc) 



Likewise the Lorentz matrix constraint (IV. 8c) is pN expanded to give 



= A 



[0](i) 



(IV. 18a) 



A 



A 



(1- 

mi 



f)-5^;i-2G(z;l-z;^) 



m2 

?'12 



GS. 



ij "'12 



G 
2^ 



3P2 



' 12 

2 — Pli + (ni2 • P2 n-12 
mi 



Gn{, 



' 12 



mi ^ 



0{c 



0(c- 



(IV. 18b) 



(IV.18c) 



These 6 constraints per body (12 in total) are comprised into a set of 12 elements — 0, A — 1,...,12, which 
enables us to construct the Dirac bracket out of them according to the standard rule, if we know the Poisson brackets 
between the various quantities which enter the constraints and the Dirac bracket, which are given by (IV.7a)-(IV.7c) 
e.g., we have the following Poisson brackets between the constrained variables: 



A 



A{k) 



A{n) 



5„,Af" 



(IV.19a) 
(IV.19b) 
(IV.19c) 



After applying all Poisson brackets, the constraints can be used in the results of the Dirac brackets. The constraints 
= are given above by (IV. 17c) for A — 1,2,3 for object 1 and A — 7,8,9 for object 2 and by (IV. 18c) for 
A = 4, 5, 6 for object 1 and A = 10, 11, 12 for object 2. The matrix Cab from the definition (IV. 9) is decomposed in 
a Minkowski part C(o)ab and a curvature correction part linear in G, C(^i)ab-, approximated by the needed pN order: 



Cai 



G, 



{Q)AB 



■ G(i 



(1)AB 



(IV.20a) 
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C 



AB 



r<AB 

<-(o) 



riAC ri 



r-iDB 

(l)CD<-^(0) 



(IV. 20b) 



Referring to [21] the Minkowski parts C(o)as and can be obtained analytically exact with the formal definition 
P? = \/ m'j + yielding: 



iQ)AB 



riAB 



(p?) 

nil -p-1 




'if 



' H.7 







mi ' 




_ Plipl] 

,0\2 






I 



m2 -p 



lij 



-ViikVijiSi 



■,-1 



mi) 



2ij 








7^ 



P2 -p-1 

"12 2ij 
Pup 1 3 

mi 














m2 ' 2U 
^ C'(^)(0 

2ikr2jl02 



(IV.21a) 



(IV.21b) 



(IV.21c) 



The notation is adapted to the one used in [21], whence the same calculation was done in the Minkowski case which 
should be included here as a limiting case. As a matter of fact we will only need the Minkowski part C(o)yis for our 
calculation, because all terms up to the next-to-leading pN order of the Dirac bracket are produced by the vectors 
to the left and to the right of the inverse matrix C^^ , which can be checked for all possible Dirac brackets, so no 
curvature terms of C^^ will contribute to the Dirac bracket to NLO. Interestingly the object Vuj also appears in 
surface terms of the stress-energy tensor algebra in Minkowski spacetime, see Eq. ( A7) in [25] . 



B. Transition to Newton- Wigner variables 

The Dirac bracket reads in reduced approximation to linear order in G (indicated by ~) 

{Ql, Q2}d ^ {Ql, Q2} - {Qu'^A}Cfof{^B,Q2} ■ (IV.22) 

A list of all possible combinations of quantities which may enter the Dirac bracket and their results is given in the 
appendix A. Assuming that we have calculated the Dirac brackets between all the variables entering the Hamiltonian 
the transition to Newton- Wigner variables Zi, pi, Si, Z2, P2, and S2, can be performed. The Dirac brackets of these 
new variables should be standard canonical fulfilling their natural standard commutation relations 

{z},pj,}u^S)dij, (IV.23a) 
I^WO-)^ = S.iS'j'^^"^ - + <S.fc5p(" - SuS'j'^^"^ , (IV.23b) 

and all other possible brackets should vanish. For the sake of completion we also list the standard Poisson bracket 
relation for the standard canonical Lorentz matrix A^^^^^\ which happens to be a pure 3-dimensional rotation matrix 
fulfilling the standard Poisson bracket relation 

while all other possible brackets with A^'l'^' being zero. The Newton- Wigner variables that are assumed to exist up to 
the considered order are fixed up to canonical transformations which opens up the possibility to choose a representation 
that leaves the momenta unchanged, pj = p/. In order to find the position and spin variable transformation we make 
a general ansatz with undetermined coefficients, e.g., for the position variable the general ansatz (excluding the 
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Minkowski term which is known exactly) reads to the pN order of c 



z\ = z{ + G 



. Plfc5'l(i)(fe) P2fc'S'l(i)(fc) '^12(Pl'ni2)<S'l(i)(A:) nj2(P2-ni2)S'i(i)(fe) 

U : 1- <,2 : 1- 43 1- M 



ri2 ri2 ri2 ri2 

nj2'5'i(i)(fc) 'S'i(fe)(0'S'i(i)(0"-i2 , , Snk)(i)Si(C}(i)Pik Si(k)(i)Si{i){i)P2k 

■ S>5 H (,6 — h (,7 — h (,8- 



?'12 



Ci 



Si(k)(i)Si(t){i)ipi-ni2)n^2 , > 'S'i(A;)(;)'S'i(i)(/)(p2-ni2)ni 



^2 



+ c 



10" 



, . '5'2(fc)(0'5'i(»)(0»^2 , ^ '5'2(fc)(i)'5'i(i)(0Pifc , ^ S2(k)(i)Sn,)(i)P2k 

+ Sll To r 412 ^2 r (,13 To 



' 12 



12 



. S2(k){i)Si(i)(i){pi-ni2)n'l2 S2(k){i)Si(i)(i){p2-'n.i2)n 

+ 4,14 r 4l5 ^2 



' 12 
12 



' 12 



12 



+ {1^2) (Cl6,...,C30) 



+ 0(c-5). 



(IV.24) 



It is worth mentioning that to first order 'covariant' and Newton- Wigner variables agree, so that the hat on the 
variables of the G-terms can also be thought of as erased when trying to cancel similar terms in the Dirac bracket 
relations. The ansatz (IV.24) contains 30 coefficients to be determined which in general fulfill no symmetries among 
themselves and are independent of each other, although most of them will be set to zero. Notice the coefficients also 
depend on mass factors and are therefore not dimensionless. In actuality it is more practical to choose a shorter 
ansatz which focuses on terms which are present in the Dirac brackets, because only those terms have to be cancelled 
in order to arrive at canonical Dirac brackets, all other terms in the general ansatz are a priori zero. The ansatz 
(IV.24) is therefore to be inserted into the Dirac brackets which contain the position variable, e.g., (A. la). It turns 
out that the fulfillment of the crucial Poisson/Dirac bracket relation (IV. 23a) is already enough to uniquely fix all 
coefficients in (IV.24), all other Dirac bracket relations that would make z canonical are then automatically fulfilled 
and serve merely as consistency checks. The coefficients for (IV.24) are then given as 
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(IV.25) 



and all other coefficients are zero. The complete transformation formula (including the pN expanded Minkowski term) 
reads 
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For the spin variable we make a similar ansatz, which we insert into (A.lj) and demand fulfillment of (IV. 23b) which 
again uniquely determines the spin transformation reading 
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(IV.27) 



The transformation of the Lorentz matrix is derived in the following section. 



V. REDUCTION VIA AN ACTION PRINCIPLE 



Spin in relativity can also be treated by an action principle [61-63], see also [6, 21, 23, 24, 64] and appendix A of 
[65]. It is indeed possible to derive the transformation to reduced canonical variables using an action approach. Let 
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us start with the effective action for two interacting spinning compact objects in curved space 

5eff = J dtLeff = J dt(^puil+P2^i2-lSlatnf-^S2abnf-H,s{zI,pJ,SKab,At'')^ ■ (V.l) 

The variables in the Hamiltonian Hcs (z/, pj, Sxab, ^f,"') are all independent from each other, no SSC is imposed yet 
on this stage. The variables in this action span therefore a too large phase space, because of the redundant 5*7(0) (i) 
and AM-^ de grees of freedom, which makes the phase space still unphysical, but will give us some insight into the 
non-reduced Poisson brackets between all the variables of the action. As usual the spin tensor Sab is defined locally 
through a projection onto a local Lorentz basis by a vierbein which fulfills the condition 

(i/e.b''9i^y ^^ab, ryab = diag(l, -1,-1,-1), ^^e^ ^r^ab ^ g,,u (V.2) 

with the Lorentz indices (a, 6, c, d) and spacetime coordinate indices (^, ..) 

(a,6,c,d) e {(0),(i)} and i€{l,2,3}, (m, z^, cr, . . .) G {0, 1, 2, 3} . (V.3) 

We need still another index label appearing in the definition of the angular velocity tensor VLab^ whose generalized 
momentum is the spin tensor Sabi which involves the Lorentz matrix Ka^ and therefore rotations and boosts with 
capital letters labeling body-fixed Lorentz indices in the sense that 

Aa^A-^, A^,A^, = 77,, with e {[0],H}. (V.4) 

Notice that KAb is time dependent. Now the definition is (the dot marking total time derivative) 

= Ka"^"^^ making f]"^ = -17''" = f^I'^^l . (V.5) 

The minus sign in front of the spin kinematic term ^Siab^f' in the action (V.l) is due to the sign convention for 
the antisymmetric angular velocity tensor fi"'' and the used signature. Indeed, variation of the action gives (dropping 
boundary terms) 



5Scs = dt 



opii J \ oz\ 

(1 o 2) 

where the variation of A"^*" was written in terms of the antisymmetric symbol 59°''' — A^"(5A"^^ and we used the 
relation 6^1°''' = jf SO"-'' + 2Qc^°'69''^'^ . The variation of the Hamiltonian with respect to the independent degrees of 
freedom of the Lorentz matrices is also achieved by usage of 69°'' according to 

dHeS Ab _ '^^cS „ABx\ b _ dHcS a \Bcx\ b _ dHcff ^ , 

°*~9A^ ~dA^^ --qJ^^ '^^B -gj^^ c]00 (V.7) 

The equations of motion are therefore given by 

= ^^''^ = 45,.[.,,],^-2^.[.A,„^, (V.9) 

and from the definition 

z\ = {z\, -ffcff} , Pli = {Pli, -ffcff} , Slab — {Slab, -ffcff} 

we can thus read off the Poisson brackets (see also (IV. 7)) 

{z},pj,} = 6,,6ij, (V.IO) 

{Slab, Sled} = SicaVbd ~ SjdaVbc + SldbVac — SjcbVad , (V.H) 
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{Slab, Afj = VbcAfa - riacAf, , (V.12) 
and all other brackets are zero. The goal is to transform (V.l) into the canonical form 

■J -J \ (V.13) 



with S'(i)(j) — eijkS(k), ^ijk — andfi^*)^^^ = Aj^.j*'' A^'^l'^^-' . Notice that with the used conventions the 

formula for the angular velocity vector fi*^*^ = —^(i) = —^eyfef^'"''"-'^'' involves a minus sign. The hat labels functions 

depending on canonical position variables {z}, a'j ^"'■'} with their generalized momenta {pn, in reduced phase 

space, meaning an appropriate SSC is imposed to get rid of 5'/(o)(i) and leading to canonical conjugate variables at 
the same time. Variation of the action is completely analogous to above calculation. Only the 4-indices a,b, . . . have 
to be replaced by 3-indices (i), (j), . . . , so the Poisson brackets (V. 10)- (V.12) translate into 

{z},pj,} = 6}5ij, (V.14) 
{Af^'\ ^|™)(")} = -,5,™AW("' + ,5,„aW(") , (V.15) 

all other zero. We used Tj(^i-^(^j) = —6ij and the antisymmetry of the spin tensor. So let us make the reduction in phase 
space explicit. We impose the covariant SSC, or to put it more exact the Mathisson-Pirani SSC with the 4-velocity 
coupled to the 4-dimensional spin tensor. As already mentioned, the reason is that the considered potentials are only 
valid for this SSC. This reduced action is then ready to be transformed to (V.13) while emerging with the Newton- 



and space parts by using the supplementary conditions fixing the frame of reference (see also (IV. 8)) 



Wigner SSC. First examine the term ^Sab^""^ (with suppressed particle label) and make a decomposition into time 



Sabu'' = ^ ^jfeW^^^O ^ '5(,)(o) =0, (V.17a) 
A['>ua = ^ aH(o)+AM(^)^ =0, (V.17b) 

M(0) 

A[o]„ = ^%. (V.17C) 

Here and in the following u"^ should be understood as given in terms of the canonical momentum p (not in terms of 
v), the pN approximate relations are given by Eqs. (IV. 16a) and (IV. 16c). Notice that we could have also chosen 

(II. 4) as SSC with p^^ already eliminated by the mass-shell constraint p^^ = ^Im^ -|-pRpPs. But again due to the 



WW 

difference between Dixon's momentum and the canonical one in the Legendre transformation it is easier to work with 
the 4-velocity in the SSCs. As the full derivation of the variable transformation formulae are quite cumbersome we 
have put the details in the Appendix B and shall present here only the key steps and results. We start with the 
insertion of the constraints (V.17a)-(V.17c) leading us the following 'naive' reduced expression of the spin coupling 
term in the action 

= -\'umi^^^''''' - ^%)wSa,"'^aW(o) + \s,,,,,^mi^) (V.18) 

with iie^)^') = Aj.j^'^^Al^K'). Notice the formal difference to the definition of VL^^'^^^^ from (^.5), Df-^'^^^^ is therefore 

not necessarily antisymmetric, which is actually an unwanted feature. After insertion of AW^'^-' by using (V.17b) and 
further algebraic manipulation we end up with the expression: 

1 „b / u'^^^u^j) u^'"^H€) \ ^''^^^'^ 1 uC^)^'^) 

-Sab^ = + S{i){k}—(7r) 'S'(i)(fc)-77vj 7^ 1- o'S'(i)(fc)-77yj ■ (V.19) 

Next thing to do is to redefine variables so that the canonical structure of (V.13) is produced. Obviously one should 
start by shifting to i^^^^'-'^ which should be antisymmetric in order to be the correct velocity variable belonging 
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to the spin tensor. This is achieved by a redefinition of the Lorentz matrix according to 
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(V.20) 



After a further redefinition of the spin tensor to the canonical (hatted) one according to 
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o, ■ ■■ 



u(u + U(o)) 



u(u + U(o)) ' 



(V.21) 



we arrive at the following reduced expression for the spin coupling term with one term, the 2^-term, left to be cancelled 
by a proper position variable shift, because this term includes a local accelaration 



iSabr!'"' = Z with z = s, 
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(V.22) 



Notice that we could have also used the principle of general covariance to arrive at (V.22), because all the equations 
up to (V.22) look the same in the special relativistic case for global Minkowski spacetime, where the round brackets 
around the local indices of the corresponding variables are erased to yield coordinate indices. The principle of general 
covariance in this case would be to rewrite the round brackets around the coordinate indices to arrive at valid expression 
in curved spacetime of general relativity. Indeed we recover the result of Hanson and Regge [21] for the transformation 
to the Newton- Wigner spin variable S and Lorentz matrix AI'l^-'^ in the special relativistic case, when the 4-velocity 
u'^ in our formulae is replaced by Dixon's momentum (II. 5). The same will be true for the transformation to the 
Newton- Wigner position variable in the special relativistic case. 

To solve for the momenta in (V.22) one has to insert the vierbein which is perturbatively calculated to the needed pN 
order, see Eqs. (IV.15a)-(IV.15c). First we make an expansion of Z in powers of (in the sense of a post-Newtonian 
approximation) 
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The goal is to find the shift of the position variable to its canonical one only approximately to linear order in G and 
to leading order in spin-orbit, spin(l)-spin(2) and spin(l)-spin(l) interaction. We insert (IV. 16a) into (V.23) and pN 
expand the result up to the order leading to 
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Next we insert (IV. 16c) into this equation yielding the approximate expression 
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We eliminate the time derivative of uf^ by shifting it on-shell (i.e. we neglect total time derivatives symbolized by 
») onto the momentum leaving us also with a time derivative of the canonical spin, which we will have to deal with 
later when we reconsider the spin redefinition. So 
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Now we are ready to return to the action (V.l), wherein we insert Eqs. (V.22) and (V.25) leading to the expression 
(for particle 1) 
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(V.27) 



3Gp2i Gni2(p2-ni2) 2Gm2n'l2Sl 



k Qii)ik) 



ri2 ri2 
This enables us to read off the position coordinate shift 
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0(c-6). 



This formula coincides with Eq. (IV. 26) when spin and position variables on the right hand side of Eq. (V.28) are 
provided with a hat in the highest pN terms (meaning all the linear in G terms here) , which is allowed when working in 
a perturbative scheme. Again the Minkowski term here is shown for pedagogical reasons only, because the Minkowski 
case can be treated exactly and in order to arrive at the pN order one has to include one higher Minkowski term in 
(V.28), which changes the coefficient when transforming it perturbatively to the left of Eq. (V.28). For the Minkowski 
case we state that we are always able to write 
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It follows from Legendre transformation of the point particle Lagrangian (u^*^ = in SRT) 
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or from the 3-components of Dixon's momentum (II. 5) which happens to be the same as the canonical one but only 
in the Minkowski case. Then the addition to the action reads 
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This whole contribution can be absorbed by redefining the position variable as 
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This is exactly the formula (B.ll) from [21]. The expanded expression yields 
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which thus matches the Minkowski terms in Eq. (IV. 26). 
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As already mentioned the spin and the Lorentz matrix need another redefinition in order to cancel the term 

~ 2m7'^i(«)0')^'i»"i"'^ from the action. This is achieved by an infinitesimal rotation o;^'-"--''' = — oi'^^'f'^ of the local basis 
so that the canonical spin and Lorentz matrices are corotated according to 
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As described in the Appendix you can read off 



, (m) , (i) 

2 mi 2 mi 



and use Eq. (V.21) to determine the final spin redefinition to our approximation: 
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Next we insert Eq. (IV. 16a) for U(o) in the second term and Eq. (IV. 16c) for u^*^ the third term of (V.36) and make 
a pN expansion up to the order c^^: 
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We are left with an insertion of the remaining 4- velocities. Again utilizing Eq. (IV. 16c) a further expansion of (V.37) 
yields 
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The last step involves providing all spin and position variables with a hat in the highest pN terms of (V.38), so that 
we end up with the transformation formula 



^PikPij Gm2PikPi] , 3Gpi-,p2fe , Gpijn'l2{P2 ■ ii-12) 



2 ml 

Gm2n'i25'i'''''Vij Gn\^S)^ 



m\ri2 



Pij 



2 mifi2 2miri2 
-(zoj) + 0(c-'). 



(V.39) 



Notice that Eq. (V.39) (or rather Eq. (IV. 27), because one has to extend the Minkowski term to arrive at the order 
c~^) is already high enough in pN order to be used for transforming effective NNLO Hamiltonians to canonical ones, 
whereas for that purpose Eq. (V.28) needs to be extended to the order c~^°, which in turn means to calculate the 
vierbein components (IV.15a)-(IV.15c) to higher pN orders. In order to get the Minkowski case right one starts with 
of Eq. (V.21) and replaces the 4- velocities by Eq. (V.29). The result can be expanded to match the Minkowski terms 
in Eq. (IV.27). 

The most important relation derived in the present section is equation (V.22). To arrive at this expression the 
transformations of Lorentz matrix and spin tensor shown in (V.20) and (V.21) were used. Notice that no approximation 
was used to arrive at (V.22), only the validity of the supplementary conditions (V.17a) and (V.17b) is required. The 
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calculation following (V.22) is devoted to the absorption of the term Z by various further variable transformations, 
most notably by a redefinition of the position. This is necessary due to the presence of the acceleration u^'^ in Z. 
However, within a perturbative context such terms in the action can be treated in a simpler way by the method 
developed in [58, 66], which in most cases amounts to an insertion of lower order equations of motion. This actually 
corresponds to an implicit redefinition of variables, so a comparison with the Dirac bracket approach is more difficult 
in this case and therefore we have not proceeded in this way here (e.g., it is likely that the canonical momentum is 
implicitly redefined). However, for an application at even higher pN orders one can take —Z as an addition to the 
Hamiltonian, after -it^*) therein was eliminated using the equations of motion. Also a mixed approach may be useful. 
One can always write 

= ^ + v% (V.40) 
TO 

where V* includes all pN corrections to this relation, see (IV. 16c) for the NLO case. Z then reads 

Z=^^lM^(Pl+v] , (V.41) 



u{u + U(o)) \m 

The first term can now be absorbed by a transformation of the position of the form 

,^^,^+ ^ . (V.42) 

mu[u + U(o)) 

while the second term is considered as a contribution to the Hamiltonian. This mixed approach thus consists of an 
explicit redefinition of the position implementing the leading order flat space transformation to the Newton Wigner 
position, followed by implicit variable redefinitions due to the insertion of equations of motion. 

The discussion in the last paragraph can even serve to modify the Feynman rules of the EFT formalism to use 
reduced canonical variables from the very beginning. First let us explain why such a modification of the Feynman rules 
provides an improvement. It is most desirable to formulate the Feynman rules in terms of reduced spin variables (either 
covariant or canonical). First, this allows for intermediate simplifications and leads to more compact expressions for 
the potentials. Second, further field variables would appear at a later stage otherwise, after the potential modes 
where supposed to be integrated out. This is due to the fact that the covariant SSC must be handled at some stage, 
which introduces the velocity in the local frame and thus further field variables, see (IV. 14); at lower orders, this 
issue is just semantics, but at higher orders it is of practical relevance. To achieve this goal one has to discuss the 
kinematic term —^Sab^"^^ in the manifestly covariant Lagrangian, which is disregarded in the Routhian approach 
as it contains no interactions with the field. ^ However, when the covariant SSC is eliminated at the level of the 
action this kinematic term turns into (V.19), which corresponds to the complicated kinematics described by the Dirac 
bracket for the covariant SSC and produces further interactions with the gravitational field via the velocity in the 
local frame. It should be stressed that due to the use of a Lagrangian all equations of motion can be obtained by 
a variational principle at any stage, without the need to resort to nonreduced Poisson brackets as in the Routhian 
approach. Still we think one should transform the reduced covariant spin to a reduced canonical one, as the kinematic 
terms (V.19) then simplify to (V.22). Incidentally this implies that the spin kinematics can be described by standard 
reduced Poisson brackets, with some advantages already mentioned in the Introduction. This procedure is essentially 
a straightforward adaption of the approach in [23] to the EFT formalism, where both the reduction of the covariant 
spin and the transition to a reduced canonical spin succeeded at the action level in a consistent and transparent 
manner (see also [24] for a more detailed exposition). 

Having this said, we propose the following steps to improve the Feynman rules. First a manifestly covariant 
Lagrangian is obtained by adding —^Sab^"''' given by (V.22) to the initial Routhian defined by Eq. (7) in [8]. Next 
the covariant SSC is inserted in the usual interaction terms of the Routhian and the spin is transformed to the 
reduced canonical one using (V.21). Finally one can use a mixed approach to eliminate the acceleration in Z by 
writing ii^^^ = e^^\u^ + e^^\u'' + e'^-'^o when — 1. Again the first term can be absorbed by a corresponding 
redefinition of the position variable and the other terms simply provide further interaction terms which must be 
included in the Feynman rules. It is left as a future task to work this out in detail. The appearance of higher order 
time derivatives poses no difficulties at this stage, as one is working with a Lagrangian. After the Feynman rules were 
applied, one may even obtain a fully reduced Hamiltonian in the usual way, if desired (i.e, by eliminating accelerations 
and higher order time derivatives followed by a Legendre transformation in the velocities). 



^ Notice that the Routhian in [8] is therefore not manifestly covariant. By adding —■^Sab^'^^ to it one can get back to the manifestly 
covariant Lagrangian contained in Eq. (1) of [8], see also Sect. Ill in [13] and Sect. 5.2.2 in [24]. 
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VI. FINAL COMPARISONS OF POTENTIALS WITH HAMILTONIANS 



In this last section we make use of the transformation formulae we have found throughout the previous sections, 
which enable us to transform all the non-reduced effective Hamiltonians from Sect. Ill to reduced ones depending on 
standard canonical variables. We are especially interested in the effective NLO Hamiltonians, which shall be compared 
with their ADM counterparts that were calculated directly within the ADM approach, see [25, 39, 42-44]. As the 
variable transformation of the covariant variables to Newton- Wigner ones has to be done in all Hamiltonians up to 
2pN order we will again get NLO correction terms stemming from all the subleading order Hamiltonians starting 
with the Newtonian Hamiltonian from (1.3). Therein we replace the position variable z by z utilizing Eq. (IV. 26) or 
(V.28). The leading-order spin-orbit correction term emerging from this procedure is labeled as ^^loso- upper 
index refers to the Hamiltonian where the variable replacement is made and lower one is reference to the pN order 
and the specific Hamiltonian which is to be corrected by that. Likewise we label all the NLO correction terms and 
get 
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The same replacements of position variables must be performed in the EIH potential giving rise to the following 
correction terms to the NLO spin-orbit Hamiltonian: 
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Next we replace variables in the LO spin-orbit Hamiltonian from Eq. (III. 6), where we first have to insert the covariant 
SSC from Eq. (IV. 17c) to arrive at a consistent expression. The spin replacement by Eq. (IV. 27) yields the correction 
term 
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Equally replacing the position variables by Eq. (IV. 26) in this Hamiltonian gives further correction terms for the NLO 
S1S2 and NLO 5^ Hamiltonians 
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whilst the replacement of the position variable in these two Hamiltonians leads to no new NLO correction terms. 
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A. The NLO spin-orbit Hamiltonian of Levi 

We adopt Eq. (III. 23) and make the transition to a canonical Hamiltonian, which depends solely on the hatted 
variables (x/, pj, S^^ ). First thing to do is to eliminate the S'l"'"-*'' variables by imposing the SSC from Eq. (IV. 17c). 
The result is the effective Hamiltonian still depending on (x/,pj, Sk)- 
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To make the transition of this Hamiltonian complete we put a hat on its variables and add up to it all the needed 
correction terms via 
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which results in 
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which is to be compared with the NLO spin-orbit ADM canonical Hamiltonian first derived by Damour, Jaranowski 
and Schafer [39], which reads 
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If both Hamiltonians are correct and are therefore to generate the same equations of motion, the difference between 
these two should equal an infinitesimal canonical transformation, which in turn involves a generator function g that 
is to be chosen appropriately as outlined in e.g., [13]. So it should hold 
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where for the generator function one can make the following general ansatz to canonically transform a NLO spin-orbit 
Hamiltonian: 
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It turns out Eq. (VI.15) can be fulfilled by choosing the coefficients to be 
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We kindly note that the same agreement of (VI.ll) with (VI.14) was already achieved by Levi [13]. However, the 
transformation to canonical variables was found in [13] by comparing with the Hamiltonian in [39], whereas here 
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we derived it from general principles. For this reason our Eq. (IV. 26) also contains a term which is irrelevant for 
translation invariant quantities like the Hamiltonian (but would be needed, e.g., for the center of mass), whereas such 
a term was omitted in Eq. (121) of [13] as it is not needed there. Notice in [39] there was also a variable transformation 
formula determined to achieve a comparison with a result which was obtained in harmonic coordinates and with a spin 
whose lenght is non-conserved. Their transformation formula (6.11), which was uniquely determined and comprises 
non-canonical and canonical transformations, contains a similar irrelevant term when comparing translation invariant 
quantities. By our derivation of Eq. (IV. 26) from general principles we also achieved direct justification of their 
transformation formula, that was left as future work by the authors. 



B. The NLO spin-orbit Hamiltonian of Porto 



We also want to compare the result of Porto with the ADM NLO spin-orbit Hamiltonian from Eq. (VI. 14). We start 
out from the non-reduced effective Hamiltonian from Eq. (III. 25) and insert the covariant SSC from Eq. (IV. 17c), 
which results in the reduced effective NLO Hamiltonian in 'covariant' variables 
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The transition to the canonical Hamiltonian is made as in the case of Levi's effective Hamiltonian by adding up the 
missing pieces and putting a hat on the variables of Porto's Hamiltonian: 
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The result is 
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Again for the difference AH^Yjjso there should exist a generator function Eq. (VI. 16) with determined coefficients. 
The coefficients for this case read 
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which means we have achieved 'on-shell' agreement with the ADM Hamiltonian. 
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C. The NLO spin(l)-spin(2) Hamiltonian of Porto and Rothstein 

For reason of completeness we also compare the NLO spin(l)-spin(2) Hamiltonian which results from the corre- 
sponding effective potential that Porto and Rothstein have calculated, with the corresponding ADM Hamiltonian. 
We refer to Eq. (III. 27) and insert the SSC from Eq. (IV. 17c) to arrive at the reduced effective Hamiltonian, which 
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reads 
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We transform this 'covariant' Hamiltonian to the canonical one by adding up to it the corresponding correction terms 
which follow from the variable transformation and replace the 'covariant' variables with Newton- Wigner ones in the 

original Hamiltonian 
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the result being 
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This Hamiltonian shall be compared with the corresponding ADM Hamiltonian which was calculated by us in [25, 40]. 
It reads 
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To generate the difference AiJ^^b S1S2 ~ ^nlo'sTs] ~ ^^^'s^rcd^ '-'^ ADM Hamiltonian to the one resulting from 
an effective potential we make another ansatz for the generator function g: 
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The coefficients can be uniquely determined and yield 
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An alternative potential in the covariant SSC was derived by Levi [9] and the procedure to obtain a fully reduced 
Hamiltonian used in the present paper (and briefly presented before in [44]) was already applied therein. Therefore 
we will not repeat it here. However, the transformation to canonical variables was found in [9] by comparing with the 
Hamiltonian in [40], whereas here we derived it from general principles. 

Notice that in [40] it was found that the result in [79] was incomplete, however, once completed it also fully agrees 
with the corresponding NLO spin(l)-spin(2) ADM Hamiltonian [8]. Though the result in [79] was derived within 
the EFT approach, it was not based on the Routhian method but on a direct insertion of a Newton- Wigner SSC.^ 
The alternative NLO spin(l)-spin(2) potential derived by Levi even agrees exactly with the corresponding ADM 
Hamiltonian if a suitable Newton- Wigner SSC is used at the action level [9]. More work is needed to understand why 
such a direct insertion of a Newton- Wigner SSC leads to correct fully reduced Hamiltonians in some cases. 



^ However, in [79] this procedure is not justified. For a proper way to directly implement a Newton- Wigner SSC in an action see [22], in 
particular Eqs. (3.7) and (3.12) therein. 
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D. The NLO spin(l)-spin(l) Hamiltonian of Porto and Rothstein 

The last Hamiltonian we want to verify is the effective NLO spin(l)-spin(l) Hamiltonian from Eq. (HI. 32). From 
it we eliminate by Eq. (IV. 17c) resulting in 
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We proceed to the canonical Hamiltonian by collecting all matching correction terms and by replacing covariant 
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variables with Newton- Wigner ones in the above Hamiltonian 
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The comparison with the ADM NLO spin(l)-spin(l) Hamiltonian was already briefly addressed [44], and we also refer 
to [44] for details on the derivation of the ADM Hamiltonian. Notice the two misprints in Eq. (32) of [44] if one 
compares it with (IV. 27). The round bracket term should be corrected to read 
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and the third term has a wrong particle label, rather it should read 
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It turned out that with the following general ansatz for the generator function 

5nlo s? = 7i(Pi-ni2)S?+72(P2-ni2)S?+73(Si-ni2)(Si-pi) 

(VI.33) 

74(Si-ni2)(Si-p2) + 75(pi-ni2)(Si-ni2)^ + 76(P2-ni2)(Si-ni2)' 

of an infinitesimal canonical transformation which belongs to this Hamiltonian, the difference to the ADM Hamiltonian 

Ai^NLo sf - {^^N, g^To sf } (VI-34) 

can be generated by choosing the coefficients to beat lower orders, this issue is just semantics, but at higher orders it 
is of practical relevance 

71 = - ^ + C'qi , 72 = , 73 = i , 74 = , 75 = , 76 = . (VI. 35) 

We have thus shown rather clearly how one can transform a non-reduced effective potential at NLO to a canonical 
Hamiltonian in a very replicable and systematically way without losing oneself in subtleties. 

VII. CONCLUSIONS AND OUTLOOK 

In this paper we have presented two different methods to transform non-reduced effective potentials, that can 
be calculated within the EFT approach, for spinning compact binary systems in general relativity to fully reduced 
canonical Hamiltonians, which then depend only on the physical degrees of freedom. The main subject we focused on 
were potentials of constrained systems with redundant spin degrees of freedom that had to be eliminated. The key 
difficulty was to achieve this scheme on curved spacetime at least perturbatively up to next-to-leading order. The first 
method involved working with Dirac brackets in Sect. IV to reduce the phase-space variables and the second method 
was to redefine variables in such a way that the action in Sect. V becomes fully reduced. Both methods yielded the 
same transformation formulae for spin (IV. 27, V. 39) and position variable (IV.26,V.28) although the action approach 
is much more transparent and provides quite general transformation formulae that possess validity in full general 
relativity given that the vierbein field of curved spacetime is known. Furthermore we could determine how the 
Lorentz matrices should transform in a gravitational field in order to arrive at a standard canonical representation, 
see Eq. (V.20). Thus our formalism is also valid for effective potentials/Hamiltonians that depend on the Lorentz 
matrices, which would describe an asymmetric behavior of the physical system in consideration. 

Interestingly the method with the Dirac brackets needs no transformation of the Lorentz matrices, which was 
expected, because the Hamiltonians are all independent of the Euler angles and the Dirac brackets operate on the 
level of the equations of motion, which are in turn generated by these Hamiltonians. In contrast the action approach 
needs the redefinition of the Lorentz matrices via Eq. (V.20) because the action depends on them while coupling to 
the spin. Due to the knowledge of the vierbein field to next-to-leading pN order we could explicitly calculate the 
transformation formulae for transforming all next-to-leading order effective potentials known to date for two self- 
gravitating spinning compact objects to canonical Hamiltonians while eliminating the spin supplementary condition 
and performing a Legendre transformation. In the near future we plan to go to next-to-next-to-leading pN order in 
the transformation formulae in order to compare the non-reduced NNLO spin(l)-spin(2) effective potential by Levi 
[15] with the fully reduced NNLO spin(l)-spin(2) canonical ADM Hamiltonian calculated by Hartung and Steinhoff 
in [48] . Clearly for this purpose it is much easier to start from the generally valid formulae found by transforming the 
action directly instead of elaborating more on Dirac bracket calculations. The only piece missing is the vierbein field 
to next-to-next-to- leading pN order, which can be calculated, e.g., from the metric in harmonic coordinates [37] and 
corrections [38] , and poses no conceptual problem. It is interesting to note that now after 37 years we finally succeeded 
with the proposal made by Hanson and Regge [21] in their conclusion section to "attempt to include gravitation in 
our formalism" initiated in [8, 22, 30, 63, 64]. Furthermore at the end of Sect. V it was discussed how Eq. (V.22) can 
be used to improve the Feynman rules of the EFT formalism by formulating them in terms of reduced canonical spin 
variables. 

One future task would be to extend this formalism for dipole (spin) interaction that we developed especially for 
weak deformation effects to higher multipoles that would account for stronger rotational deformation effects which 
will play an important role when the binary system is very close to the merger phase. Right before the merger there 
is also no more guarantee for the spin length to be conserved which would also be an attractive topic to deal with, 
because so far there has been no algebraic treatment of a change in the spin length, only analytically starting with 
Teukolsky's analysis [80]. 
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Appendix A: Dirac brackets 

Here we list all possible combinations of quantities to enter the Dirac bracket and their results except the ones 
containing the A-matrices. The reason is that all Hamiltonians of interest in this work do not own any A-dependency, 
so there is no need to canonicalize the A-matrices explicitly here with the Dirac brackets besides all the other variables 
{z} ,pji, S'^), which really appear in the Hamiltonians. This is different when eliminating the SSC and transforming 
to canonical variables on the action level, see Sect. V. Then the first thing one must do is to transform the A-matrices 
to canonical ones being the hatted A-matrices. The results of the Dirac brackets read 



Zi,Pij}d ^ + — 2"i2Pife'5i 



2mf 



2mi 



ni2P2jS 



{^)ik) 



1 



2toi 
G 



(p2-ni2)5i 



2toi 



'^i2(P2-ni2)rii2'S'i 



(i)(fc) 



+ ZT-^- H"l2'^12i'l ^2(fc)(0-i'i ^2{j)(k) 

+ ( 3'^i2^l2^i ^iw(o - ^1 ^lum )+o{c ) 



S^i cO')('-')l 
1^1, "^i Id 



(fe)(0 



ml 

2pi[fcS'i(j)](/)PiH 



™ q(i)(j) 
TO^ri2 



(A.la) 



(A.lb) 



G 



35'i(j)[(fe)P2j] 



4m5 



miri2 \ mi 

'2.GSi(i)[(j)S2(k)](i)ni2i 2Gm2S'i(i)(;)S'i(/)[(j)ni2fc] 



Sl(i)[(j)Plk] + (P2-ni2)S'i(i)[(fc)ni2j] 



n' 12 



(A.lc) 



1^17 2^2 j-D — I'J'Jl + ^1 "-12fc"i2 

2miri2 V 



G 



2m2ri2 



35- 



s. 



ni2kni2 



Oic 



\z\,^2 Id 



{PUiPljSl- 

XPu-.'^x |r 



3G5'2(j)[(fc)P2j] G''5'2(;)[(fc)P2j]"-12"-12 2G5'2(j)(;)5'i(j)[(fc)ni2;] 



n^2?'i2 



2G5'2(fe)(/)5i(,)[(,)ni2jj _5 



m-Yr\^ 



"^2^'l2 



m\r\2 



{Pli,P2i}D = {P2i,P2j}D =0, 

G 



m\r\^ 



vP2-ni2)5i(i)[(j)Pn;] 



m2 



Pil]S^k)]{i) 2 piin\2 ~ 5p2in\2 + P2tn^2 - 3(p2 ■ni2)ni2n-l 



mi 



2Gpi[kSl(j)](m) ( , I (0( 
™ „3 1 «J"-12"-12'-'2 



mir^2 



s. 



2(»)(m), 



(A.ld) 

(A.le) 

(A.lf) 
(A.lg) 



(A.lh) 



33 



2Gm2Pi[kSi{j)-\(i) 



V 12 



(3'^i2'^i2'S'i(m)(/) - + C' (c ®) , 



{ph,4^'^^'^}d--{ph,^P^}d (1^2), (A.li) 
8G'TO2Pi[i'5'i)(j)][(fe)Pi/] GG 



^2^^2 ' ^^^^^ (P2[j5'i(i)]p)Pu] +Pi[jS'i(,)][(fe)P2/]) 

2G'pi[j .S'i(t)] [(fc) ni2i] (p2 ■ ni2 ) ^ 2G'ni2[j5'i(i)][(fc)Pi;] 
TOiri2 "iiri2 

4G 

+ zr-;r (pi[j'S'i(fc)][W'5'20)](m)'^S +Pl[J'S'lW][(/)'S'2(fe)](m)?^^2) (A-ij) 

"'1' 12 

+ ('^i2[;'S'i(fc)](™)Pi[jS'i(,)](„i) + ni2bS'i(,)](™)Pi[fcS'i(,)](„)) 

'''■1' 12 

+ 0(c-«) , 

{5'f''-'\S'^''^^'^}D 2 (P2[feS'2(;)](™)ni2[jS'i(j)](,„) +?l5^S'i(„)[(j)S'2(,)][(fe)P2/]) 

'"■2' 12 

4G 

+ Z— 3- (Plb'5lW](m)»^12[i'52(fc)](m) + Pi [^ '5'l(i)] [(;) '52(fc)] (m) ^2) 
'"-1' 12 

+ O (c-*^) . (A.Ik) 
All remaining Dirac brackets follow via exchange of particle indices 1 with 2 in all terms. 

Appendix B: Full derivation of the transition to canonical variables via an action principle 

In this section we shall elaborate in great detail on the reduction of the action and the resulting transition formulae 
to canonical variables presented in Sect. V. The first step involves the insertion of the constraints (V.17a)-(V.17c) 
into the spin coupling term in the action, which yields 
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Because of antisymmetry of the spin tensor the first term with the 4-velocities is irrelevant for (B.3). We are left with 

Next thing to do is to redefine variables so that the canonical structure of (V.13) is produced. Obviously one should 
start by shifting to which should be antisymmetric in order to be the correct velocity variable belonging 

to the spin tensor. To find the correct redefinition it is useful to make an ansatz for the intrinsic redefinition of Lorentz 
matrices with an unknown function ^ according to 
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transforming A^'l^^'' into a 3-dimensional rotation matrix yielding 
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Taking the limit = the function ^ should still be regular because this state corresponds to the rest frame limit 
which is therefore an adequate physical limit to take. This argument leaves us with the sole solution 
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So the Lorentz matrices undergo a redefinition 



35 



Next we take the time derivative of this expression 
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The Symbol [=] means we neglect all terms symmetric by interchanging i j and keep those being antisymmetric 
by this interchange, because only those will contribute when projected on to the spin tensor <S'(i)(j) or equally an- 
tisymmetric expressions. We insert above expression into Eq. (B.5), which leads to a rather long expression. We 
define 
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and expand it 
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Hence the spin coupling term is reduced to the expression 
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To solve for the momenta in (B.26) one has to insert the vierbein which is perturbatively calculated to the needed pN 
order, see Eqs. (IV.15a)-(IV.15c). First we make an expansion of Z in powers of (in the sense of a post-Newtonian 
approximation) 
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Next we insert (IV. 16a) into (B.27) and pN expand the result up to the order leading to 
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After that we insert (IV. 16c) into this equation yielding the approximate expression 
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We eliminate the time derivative of w^"''' by shifting it on-shell (i.e. we neglect total time derivatives symbolized by 
«) onto the momentum leaving us also with a time derivative of the canonical spin, which we will have to deal with 
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later when we reconsider the spin redefinition. 
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Now we are ready to return to the action (V.l), wherein we insert Eqs. (B.26) and (B.29) leading to the expression 
(for particle 1) 
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This enables us to read off the position coordinate shift (for the Minkowski case only the leading order term is shown) 
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+ 0(c-6). 



1 (i) 

The spin and the Lorentz matrix need another redefinition in order to cancel the term ~ 2^ ^i{i)U)P'^'^ i from the 
action. This is achieved by an infinitesimal rotation of the local basis so that the canonical spin and 

Lorentz matrices are corotated according to 
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We identify 
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and use Eq. (B.23) to determine the final spin redefinition to our approximation: 
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Next we insert Eq. (lV.f6a) for U(-o) in the second term and Eq. (lV.f6c) for u*^*' the third term of (B.35) and make 
a pN expansion up to the order c^^: 
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We are left with an insertion of the remaining 4- velocities. Again utilizing Eq. (IV. 16c) a further expansion of (B.36) 
yields 
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The last step involves providing all spin and position variables with a hat in the highest pN terms of (B.37), so that 
we end up with the transformation formula 
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